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I ! Abstract 

^ I Some inverse problems in high-energy physics, neutron diffraction and NMR spec- 

^ I troscopy are discussed. To solve them, the Fourier integrated transformation method and 

^ • the Maximum Entropy Technique (MENT) were used. The integrated images of exper- 

Q ■ imental distributions are shown to be informative when determining the space-time pa- 

^ . rameters of a particle generation zone and when analysing blurred spectra. The efficiency 

^ ! of the above methods was checked by comparing relevant results with the results obtained 

^ I independently. 

^ ■ An inverse problem arises where the information required can not be obtained through 

ly-^ ! direct measurements. The majority of inverse problems in nuclear physics is connected to 

O I interpretation of experimental data on scattering high-energy particles on targets. The 

2 I scattering pattern is used to determine the structure of colliding particles, interaction 

T— I ' potential and the dynamics of processes occurring upon collision. In some cases, a nucleus 

^ ■ or a particle is understood as a probe with known magnetic moment, wave length and 

^ . energy level system. Its interaction with a substance and a field is used to assess the 

I ! concentration of the substance, the characteristics of the medium or the configuration 

Qj \ of the field. Such inverse tasks are typical of NMR and NGR spectroscopy, neutron 

^ ' diffraction, activating analysis, emitted tomography and other applications of nuclear 

^ ■ physics. 

^ . In this paper, some inverse problems arising in experiments in high-energy physics are 

^ ! discussed [|I], 0, |, ^, |^. To solve them, filtration methods are used. Experimental data 

" " ' are subjected to Fourier transformation to obtain a result without a substantial use of a 

priori restrictions |^ ^. This property is especially important when selecting models of 
processes in cases when the experimental distributions are compatible with several models 
and the application of approximation methods is ineffective. 

One question, which arises when solving an inverse task, is whether the results ob- 
tained are proven. This question can be answered by comparing the results obtained by 
independent methods. To illustrate such a comparison, diffusion spectra were used as 
an example, their resolution being increased instrumentally and a posteriori through the 
Maximum Entropy Technique (ENT) 0, ^ . 

The goal of our study was to show that the same type of mathematical and algorithmic 
approaches, based on integrated transformations, can be used to solve some inverse tasks 
in nuclear physics. 

SPACE-TIME CHARACTERISTICS OF THE MULTIPLE BIRTH OF PARTICLES 

Using Podgoretsky and Kopylov's approach, light can be thrown on the space-time 
characteristics of the multiple birth of particles by studying correlations of identical parti- 
cles with close 4-impulses [Q. The probability W{qo, q) of detecting two identical particles 
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with 4-impulses Pi = {si,pi} and P2 = {£2,^2}, emitted by two point sources with the 
coordinates {ti,fl} and {^2,^^}, is expressed by the formula: 

I cos[g(r-l-r-2)-go(ti-t2)] 

where q = Pi — "fi, Qo = £1 — £2, t is the hfetime of the sources. For identical particles of 
identical energy go = 0, a relationship (1) is determined by a spatial variable R = 7\ — r^: 

W{q) = l + cos{qR). (2) 

The above formulas hold for identical particles that obey Boze-Einstein statistics, e.g. the 
TT-mesons of large impulses. Let us assume that f{R) is the normalized spatial distribution 
of the sources of particles on the R. Then the formula (2) for the probability W{q) has 
the form: 

W{q) = 1 + y /(^) cos{qR)dR. (3) 

In this case f{R) = f{—R), because two identical particles are indexed by two equivalent 
methods. The function W{q) does not depend on the sign R and the integrated addend 
in (3) is the Fourier image F{q) of the function f{R) 

W{q) = l + F{q). (4) 

Thus, the experimental relationship W{q) is used to obtain the Fourier image of the 
distribution function f{R) 

F{q) = Wiq) - 1 (5) 

and to use it, in turn, to calculate the distribution function f{R) by inverse Fourier 
transformation. 

In collisions of identical particles, experimental data can directly provide information 
on the distribution density of the sources p{r) in the generation field, provided the sources 
of the particles are independent The Fourier image G{q) of the function p(r) is 
calculated using the formula 

G{q} = y/W{q} - 1 (6) 

and the function p(r) is estimated by its inverse Fourier transformation. The three- 
dimensional function p(r) can be represented for convenience as cutting of the generation 
zone by a system of parallel planes or a projection of the generation zone to a plane 
perpendicular to the preset direction. The condition p(r) > restricts experimental 
distributions to a certain class of functions. Such a restriction is automatically removed 



if the distribution of the sources is reconstructed by its Fourier image using ENT |T0[ . 

A similar technique of selecting identical particles with g = to estimate the distri- 
bution function of the sources on the variable t = ti — 12 is impossible as it results in the 
value go = 0, which follows from a condition relating the masses, energies and impulses 
of identical particles. When discussing a general case, we will not restrict the go value. 
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Let us introduce the normalized distribution function t,r) and re- write the formula 
(1), considering a trigonometrical identity: 

cos{qR — qot) — cos{qR) cos(go^) + sm{qR) sm{qot) (7) 

as 

cos(g7?) cos{qot)dRdtdT 



W{qo,q)^l + J h{R,t,r)- 



+ 



+ I HRM ■ (8) 



The second integral addend in the formula (8) is equal to zero because the methods of 
indexing of a particle in a pair are indentical, the function fi[R,t,T) is even and the 
function of a sine on the variables R and t is odd. 

In the assumption about the independence of spatial, temporal and relaxation distri- 
butions, described accordingly by the normalized functions f{R), x(i) and (p{t), expressed 
by the condition 

h{R,t,T) = fiR)-xit)-^iT) (9) 

formula (8) is re-written as 

W{qo, Q) = ^ + J fiR) cos{qR)dR J x{t) cos{qot)dt J J^^^^y- W 

The first two integrated co-multipliers in (10) represent the Fourier images F{q) and X{qo) 
of the normalized spatial f{R) and temporal x(i) distributions. The third integrated 
multiplier does not depend on the magnitude of qo, which follows from the condition of 
normalizing of the function ip{r) 

oo oo 

J ip{r)dr = J ip{qor)d{qor) = 1. (11) 



oo 

and can be written as constant C — [ '^^^^ ^ ^ > 0. The formula (10) then has the 

J l + iqory 



form 

W{qo,q)^l + CF{q)X{qo) (12) 

can be used to estimate the function CX{qo) from the experimental relations W{qo, q) 
and W{0, q) 

- — - iy(o,g)-r ^^^^ 

The distribution function of the sources x(t) is determined by inverse Fourier transfor- 
mation of the function CX^qo) and subsequent normalization by unit. 
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In this problem, the integrated images of spatial and temporal distribution of the 
sources of particles are present in the data obtained by the interference experiment, and 
can be revealed if the data are represented in a special manner. In some cases, an artificial 
transition from experimental distributions to their integrated images and nonlinear oper- 
ations with them can be useful. There is a class of such operations that can be done in 
the case of inverse transition to obtain more informative estimates of data than in initial 
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RETRIEVAL OF BARION RESONANCES IN EFFECTIVE MASS SPECTRA 

The aim of our study was to investigate two-particle correlations of the reaction prod- 
ucts of multiple birth at high energies that are observed in the distributions of pairs of 
secondary particles on the effective mass Mg// and provide information on compound sys- 
tems or resonances formed during particle generation. Of interest for the study of collision 
dynamics are A-isobars, barion resonances that form a connected short-lived system from 
a nucleor and a vr-meson. This interest is due to the fact that the birth of such resonances 
on nuclei upon collision of hadrons with the nuclei at low energies can provide the main 
channel for elastic scattering. Isobars can give information about interaction with other 
fragments of the nucleus and enable us to experimentally study the colour degrees of 
freedom. 

The purpose of the study of the effective mass spectra (vr^p) pairs in n~p and vr'^^C 
interactions at P.^- = iOGeV/ c is to determine the structure of the spectra and to relate 
the birth of barion resonances and a special case of the multiple generation of hadrons, 
such as birth of cumulative particles in 7r~^^C interactions 0]. 

The spectra of the effective masses of pairs were plotted on the basis of non-elastic 
11688 TT'p and 8642 7r~^^C interactions obtained by analysis of stereo photos taken by a 
two-metre propane camera of LHE JINR irradiated by a bundle of of vr-mesons from the 
Serpukhov accelerator |]12| . The effective mass of pairs was calculated using the formula: 



Meff = sjml + ml + 2{EpE^ - ppp^ cos 6) (14) 

using the masses of proton and vr-meson rUp and their total energies Ep and E^^ the 
absolute values of impulses Pp and p.,, and an angle 9 between impulses (in a laboratory 
system). The impulse of the protons selected varied from 140 to 700 MeV/s. 

To establish a relation between the birth of A isobars and the cumulative irradiation 
of TT-mesons in interactions, two schemes of selection of secondary vr-mesons in the event 
were used. In the first scheme, all vr-mesons of the event were considered. In the second 

scheme, vr-meson with a maximum cumulative order value in the event was excluded. The 

E — p// 

cumulative order of vr-meson was determined in accordance with the formula ric = 

rriN 

through total energy E, longitudinal impulse p// of vr-meson and the mass of the nucleor 
niN- The structure of the effective mass spectra vr^p pairs was analysed using two data 
processing methods based on Fourier transformation. In the former method, we revealed 
the periodicities of the dependence of the module of the Fourier image of a symmetric 
spectrum on angular frequency. These, in turn, were used to determine the positions of 
structural components corresponding to resonance masses 0. The results obtained by 
this method are shown in Figure 1. The positions of the specified characteristics of the 
spectra are shown by arrows, position errors by segments above the arrows and the value 
of the mass of known resonances by a number. 
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In the latter method, the idea to control the widths of the components of the effective 
mass spectrum was used p. The widths of the components in the estimation of the spec- 
trum were diminished by multiplying the Fourier image of the spectrum by an exponential 
function. Upon inverse Fourier transformation of such an image, we the resultant estimate 
had a higher resolution than the initial spectrum and exhibited some characteristics that 
could be compared in effective mass to barion resonances Q]. 

The results obtained by this method are presented in Figure 2 which shows experi- 
mental distributions — ^ '^^^^ and calculated estimates of f{Meff) coresponding to the 

pairs TT^p in the former scheme of selection of vr- mesons of the event (a and b), and ir^p 
pairs in the latter scheme of selection (c and d). The main peak of the structures detected 
corresponds to a mass of 1.232 MeV/c2 and additional peak corresponds to mass 1.650 
MeV/c2. When cumulative vr-mesons were eliminated, the estimates of the effective mass 
spectrum declined in a uniform manner by ca. 30% without changing the form. 

The characteristics obtained in the effective mass spectra of ir^p pairs were compared 
to A+o+(1.232), A+o+(1.650), A+o+(1.670), A+o+(1.910) isobars in n'p interactions and 
to A+o+(1.232), A+"+(1.650) isobars in vr^^^C interactions. Analysis of the average multi- 
plicity of 7r+ and 7r~-mesons in tt^^^C interactions at P^- = 40 GeV/c, made with regard 
for the background of the spectra has shown that elimination of any tt~ meson in the 
event decreases by about 30% in the distributions —^wjp'- "^^^^ result is believed to 
show that the generation of cumulative vr-mesons and the birth of barion resonances are 
two independent processes. 

RAPIDITY-BASED PARTICLE DISTRIBUTION STRUCTURE 

According to modern concepts of interaction between high-energy hadrons and nucle- 
ors and nuclei, secondary particles dominantly arise during the hadron transformation of 
fast quark parton into a jet of particles [|l^. It is convenient to study this jet hadron 
transformation, seen in short-range correlations of interaction products, in the rapidity- 
based distributions of reaction products, observing the grouping of particles relative to 
several characteristic rapidities. 

Rapidity is determine as y = (1/2) ln((£' + p///{E — p//)) where and p// are energy 
and the longitudinal impulse of a secondary particle. 

Rapidity-based distributions are typically plotted as smooth curves without visible 
characteristics, which makes it difficult to select jets. To detect the characteristic rapidi- 
ties of the multiple generation of vr^-mesons in Ti~p and vr^^^C reactions at P^- = 40 
GeV/c the method, used earlier to decrease the widths of distribution components, was 
applied 0. 

Input data and the results of the application of this method are shown in Figure 3. 
The rapidity-based distributions of Tr^-mesons in these reactions (a and b) are highly 
blurred, but their contrasting estimations (c and d) are structurally complex. In Ti~p 
interaction, the estimate of rapidity-based distribution of 7r+-mesons is a doublet with 
peaks y ~ 1.4 and y ~ 2.8. For 7r~-mesons, it consist of three peaks with y ^ 1.2, 
y ~ 2.8, y ~ 4.7. The estimates of rapidity-based distribution in n'^'^C interactions for 
TT """-mesons contains three peaks with rapidities y ~ 1.4, y ^ 2.8, y ~ 3.3 and three peaks 
with rapidities y ~ 0.9, y ~ 2.5, y ~ 4.4 for Tr^-mesons with a highly diffused second 
peak in the rapidity interval of 1.8-3.2. 
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The results obtained can be used to compare the characteristic rapidities of secondary 
TT-mesons over the ranges 1.0-1.6; 1.8-3.2; 4.0-4.4 to the areas of fragmentation of a 
target, a colhding vr-meson. The number and values of the characteristic rapidities of the 
process depend on the mass of the target and the charges of secondary vr-mesons. 

INCREASING THE INFORMATIVE CAPACITY OF SPECTRA BY MENT 

A solution can be increased by diffusive spectra by the MENT. To increase the solution 
of experimental spectra, a variant of ENT with trial diffusion functions is used ^ . This 
variant is similar to Fourier method for regulating distribution widths, but has the best 
solution and stability to the noise JTI . 



In this study, MENT was used to reconstruct the inner structure of diffusive neutron 
diagrams and the NMR spectra of solid state. These spectra typically provide a tool for 
recording spectra that differ in resolution. This allows to compare the results obtained 
by processing the original diffuse low-resolution spectra with those obtained by recording 
high-resolution spectra and to find out whether diffuse spectra can be processed. 

Neutron diffraction is an effective tool for the structural study of a substance due to 
various properties of a neutron and the comparability of its wave lengths and energies of 
neutrons to interatomic distances and excitation energies in solid and liquid states. 

A convenient scheme of neutron diffraction, compatible with impulse reactors of neu- 
trons, is the reflection of neutrons from an atomic plane at a constant angle with de- 
composition of reflection intensity in a spectrum on the wave length of neutrons or the 
flight time by neutrons of the defined base. When developing this technique, Fourier 
diffraction was accomplished with the help of a scheme of Fourier breakers of a beam 
on impulse reactors of neutrons. Neutron Fourier diffraction is characterized by a high 
spectral solution and by smaller level of a hum. 

Fig. 4 shows neutron diagram of powder YBa2{Cu, Fe)30Q3Yi2 registered with usual 
solution by a time of flight diffractometer (a) and neutron diagram of this sample measured 
on the IFR-2 reactor (LNF JINR Dubna) with high solution by Fourier diffractometer 
(b) p. Outcome of treatment of diagram (a) by the MENT: estimation and its touch 
diagram are gave on (c). 

The comparison of these graphs suggests the similarity of structures represented by 
the MENT touch diagram and Fourier diagram. 

NMR spectra provide important chemical, biophysical and medical information. 

Magnetic resonance is characteristic of nuclei with a spin differing from zero placed in 
a magnetic field. It is observed as the precession of the magnetic moments of the nucleus 
in a magnetic field. The macroscopic vector of magnetization can be recorded during the 
interaction of moments with electromagnetic radiation containing frequency components 
near Larmor frequency c^o or under the influence of fast changes in the orientation of an 
external magnetic field. For this purpose, in addition to a static magnetic field, a radio 
frequency field with Larmor frequency ujq is applied perpendicularly to direction of static 
field or several types of RF- impulses (called 7i/2, n, 3n/2 and 27i in accordance with 
the angle of turn of the magnetization vector in a rotating coordinate system) are used 
to initiate transition processes in a system of nuclear spins. A signal of free collapse of 
induction after the cessation of the action of tt/2 of impulse can be used, as a result of 
inverse transformation of Fourier, to obtain data on the characteristic frequencies of a 
system of nuclear spins as well as spin-lattice and spin-spin relaxation times. 
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NMR spectroscopy can conveniently be applied to liquids the spectra of which are 
represented by narrow spectral lines. The NMR spectra of solids, recorded by the forced 
precession method, are, on the contrary, diffuse and are hard to interpret. The reason 
for that is that the averaging of dipole-dipole spin-spin interaction is less efficient than 
in liquids. Artificial averaging became possible after using special sequences of 7r/2 and 
TT impulses . The NMR spectra of solids then became similar to those of liquids. 

Applying MENT to a low-resolution NMR, we can compare the estimate of the spec- 
trum obtained with a high-resolution NMR spectrum. In this example, our purpose was 
to check up the ability of the MENT algorithm for the diffuse spectra of solids that have 
a constant noise level. 

Figure 5 shows the NMR spectrum ^^C of adamantane, an organic compound (a) 
recorded by the forced precession method (curve 1) and with the help 7r/2, vr sequences 
of impulses (curve 2) Also shown in the figure is a MENT estimate (curve 3) with its 
touch diagram (curve 4) (b) obtained by processing the diffuse NMR spectrum. 

Characteristic of a processed spectrum is its arbitrary position on the ordinate axis, 
i. e. the lack of information on the constant component of hum noise. It is usually difficult 
to process and interpret such spectrum. Indeed, the touch diagram of the MENT estimate 
of this spectrum shows a series of peaks mot of which have a background origin. However, 
there are typically two intensive peaks that reflect the true structure of the high-resolution 
NMR spectrum recorded. 

This example shows that the hum noise of a diffusion spectrum should be accurately 
estimated to draw correct conclusions about its structure. It also shows that the position 
of the most intensive components of a diffusion spectrum can be determined even though 
there is no information on the constant constituent of hum noise. 



Conclusions: 

1. Based on the interference approach, a method for reconstruction of the space- 
time characteristics of a particle generation zone in collisions of high-energy particles 
is proposed. This method can be employed to determine the density functions of the 
space-time distributions of the sources of particles from the experimental distribution of 
pairs of secondary identical particles with respect to a difference in their impulses and 
energies, thus solving an inverse problem. The statistical provision of modern experiments 
of high-energy physics favours the application of this method. 

2. The application of the methods of Fourier analysis of diffusion spectra to search for 
barion resonances and characteristic rapidities in the reactions of tt~p and 7r~^^C interac- 
tions at Ptt- = 40 GeV/c has made it possible to elucidate the structural characteristics of 
experimental distributions and to compare them with known A isobars and zones of frag- 
mentation of reaction products in the rapidity space. The results obtained by analysing 
the spectra of the effective masses of pairs using various schemes of selection of 7r-mesons 
in an event show that the formation of A isobars and the emission of cumulative vr-mesons 
are independent processes. This approach can be used to search for various resonances, 
multi-quark systems, exotic mesons and purely glue states [|1^ . 

3. In inverse problems arising in applied nuclear physics in neutron diffraction and 
NMR, MENT-based methods of analysis of diffuse spectra gave results comparable to 
those obtained by high-resolution instrumental methods. A high informative capacity 
of the MENT methods was supported by processing simple and complex spectra in the 
absence of data on the constant constituent of hum noise. 
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